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The generation of nanoscale square and stripe patterns is of major technological impor- 
tance since they are compatible with industry-standard electronic circuitry. Recently, a 
blend of diblock copolymer interacting via hydrogen-bonding was shown to self-assemble in 
square arrays. Motivated by those experiments we study, using Monte Carlo simulations, 
the pattern formation in a two-dimensional binary mixture of colloidal particles interacting 
via isotropic core-corona potentials. We find a rich variety of patterns that can be grouped 
mainly in aggregates that self-assemble in regular square lattices or in alternate strips. Other 
morphologies observed include colloidal corrals that are potentially useful as surface tem- 
plating agents. This work shows the unexpected versatility of this simple model to produce 
a variety of patterns with high technological potential. 

PACS numbers: 

It has been a challenge in nanotechnology to produce highly ordered structures by controlling 
the position of nanoparticles over an extended length scale [1]; it is not easy to achieve this goal 
using direct micro- and nano-fabrication since such processes are prohibitively expensive and time- 
consuming below a certain length scale. As a result, the search of particles on the mesoscopic 
scales that self-organize into potentially useful structures by virtue of their mutual interactions is 
extremely important [2J, [3J. The ability of soft-matter systems to self-assemble in a surprisingly 
large variety of periodic arrangements has lead to use it as building blocks of bottom-up nanofabri- 
cation processes. This is a large and rapidly growing field of tremendous technological potential and 
fundamental interest that has been feeded by the continuing progress in the manipulation of the 
interaction potentials between nanoparticles ^-pj. In particular, the self-assembly of nanometer- 
length-scale patterns in two dimensions is currently of interest as a method for improving through- 
put and resolution in nanolithography [9j; however, traditional self-assembling approaches based 
on block copolymer lithography spontaneously yield nanometer-sized hexagonal structures that 
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are incompatible with the square arrangements used in industry-standard circuits [8]. That is the 
reason why among the targeted morphologies, stripes and square arrays are particularly important 
for the manufacture of microelectronic components since they would enable simplified addressabil- 
ity and circuit interconnection in integrated circuit nanotechnology. There have been theoretical 
efforts to tailor interparticle interactions to spontaneously produce target many-particle configura- 
tions for single component systems [2j , [3j but the optimized potentials found are difficult to fit by 
available interactions. Recently, an experimental technique based in a blend of diblock copolymers 
was developed to produce nanoparticles that assemble into square arrays [8j. The blend consists 
of two kinds of diblock copolymers, A-B and B'-C, being the A and C blocks mutually repulsive 
and incompatible with blocks B and B'. Block B contains small numbers of groups that form hy- 
drogen bonds to complementary groups in block B' (Fig. la). The attractive interactions between 
complementary hydrogen bonds suppress macrophase separation in favor of microphase separa- 
tion, thereby producing large-scale assembly of nanoscale features. By controlling the amount of 
hydrogen bonding units, the molecular weights and compositions of the block copolymers, diverse 
families of ordered structures are achieved, including square arrays of cylinders that are the re- 
sult of the competition between the different interactions (see Fig. 3 of Ref. [8J). Inspired by 
these experiments, we propose a simple model system that reproduces such structures as a result 
of a binary mixture of nanoparticles interacting via isotropic short ranged competing potentials. 
The model we propose treat the interactions between the block copolymers in an effective way 
that takes into consideration the entropy of the brush-like coronas surrounding each cylinder via 
a soft repulsive shoulder and the hydrogen bonding via a square well attraction. Such simple pair 
potentials are often used to describe effective interactions among substances with supramolecular 
architecture [IOj,[IlJ. We consider a two-dimensional (2D) assembly of particles consisting of a 
hard core surrounded by a soft corona which is repulsive for particles of the same species and 
attractive otherwise. The simplicity of the model allows to predict, based on geometrical argu- 
ments, under which circumstances a given family of structures is obtained. Among the rich variety 
of patterns encountered, we can highlight the formation of square arrays and stripe phases for 
equimolar mixtures, and the formation of colloidal corrals for asymmetric mixtures. These results 
suggest a strategy for producing a range of self-assembled structures which could be exploited for 
use as templating or directing agents in materials syntheses. Although we are primarily interested 
in the crystalline structures at low temperatures, we find that at higher temperatures the system 
presents a transition from an isotropic fluid phase to a disordered fluid-like stripe phase. 

Our model is not limited to the above mentioned diblock copolymer blend since core-corona 
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architectures are also present in numerous physical systems such as dendritic polymers, hyper- 
branched star polymers, etc. Among these are, for instance, colloidal particles with block- 
copolymers grafted to their surface where self-consistent field calculations lead to effective interac- 
tions that can be modeled by a square-shoulder potential [12j. Such interactions can be controlled 
by adjusting the length, species, the grafting density of the grafted polymers, the quality of the 
solvent, the density and location of the hydrogen bonding, etc. Numerical simulations have shown 
that single component softened-core repulsive potentials may give rise to strip phases [l2]-[r7j and 
periodic structures that are explained in terms of the competing interactions between the hard 
core and the soft shoulder. Our aim is to study the influence that the introduction of a second 
component has in the domain formation. The second component interacts attractively with the 
first one in order to suppress macroscopic phase separation. This attractive interaction models the 
hydrogen bond attraction between B and B' segments of the diblock copolymer system of Ref. [8j. 
When this system forms cylindrical micelles, as shown schematically in Fig. lb, then our model 
treat them as nanoparticles with a core-corona architecture where the hard core represents the 
blocks A and C and the soft corona models the effective interactions between the blocks B and 
B'. Therefore, our system consists of a binary mixture of particles interacting through a radially 
symmetric pair potential composed by an impenetrable core of diameter ctq with an adjacent square 
shoulder with range Actq (Fig. Ic), z.e., 

oc, r < ctq, 

±e, ao < r < Actq, (1) 
^0, r > Aao, 

r being the pair distance. Particles of the same component interact through the repulsive shoulder e 
whereas particles of different species interact through the attractive potential well of depth — e. The 
model intend to represent the interactions of the block copolymers of Ref. [8j as long as they form 
cylindrical micelles (Fig lb). At large distances, the particles do not overlap and the interaction 
vanishes, the repulsive shoulder models the steric repulsion between blocks B or between blocks B' 
due to the overlap of the brush-like coronas, and the attractive well models the hydrogen bonding 
between blocks B and B'. Finally, at small separations penetration of the compact cores is very 
unfavorable and gives rise to the hard-core repulsion. The simple functional form of the interaction 
potential not only capture the essential features of colloidal particles with core-corona architecture, 
it also offers many computational advantages and allows to understand using simple geometrical 
considerations the system's self-assembly strategy [18j. Standard Monte Carlo simulations based 
on the canonical ensemble (NVT simulations) in a square box of side L with periodic boundary 
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conditions have been carried out using the Metropohs algorithm. We have used ctq and e as 
length and energy units, respectively, and we have studied the pattern formation dependence on A, 
reduced temperature T* = fc^T/e, where ks is the Boltzmann's constant; reduced number density 
p* = Ncr'^/Lp' ^ and relative amount of each species x — N2/N1, with Ni the number of particles of 
species i and N = Ni-\-N2 the total number of particles. Simulations are performed with N = 1000 
particles (runs with N = 4000 were done for the largest value of A). In all cases, the system is first 
disordered at high temperature and then brought from T* = 1 (T* = 10 for the largest value of 
A) to the final temperature T* = 0.1 through an accurate annealing procedure with steps of 0.01 
(0.1 for the largest value of A). An equilibration cycle consisted, for each temperature, of at least 
2 X 10^ MC steps, each one representing one trial displacement of each particle, on average. At 
every simulation step a particle is picked at random and given a uniform random trial displacement 
within a radius of O.Sctq. 

A variety of interesting structural features of this simple model at low temperatures are in 
evidence in Fig. 2, where a few of them are exhibited. Panel (a) shows a representative example of 
the spatial configuration that the system adopts for A = 1.5, p* = 0.5, and x = 1. We observe that 
the system self-assemblies forming an aggregate made of an alternating square array of particles 
with a lattice parameter determined by the range of the soft coronas (inset Fig. la). The system 
adopts this array in order to maximize the number of favorable overlaps between particles of 
different species avoiding at the same time unfavorable overlaps between particles of the same 
type. Thus, this simple model reproduces qualitatively well the main characteristic of the block 
copolymer system of Ref. [8j, i. e., the fundamentally and technologically important result that 
this system self-assemblies forming square arrays. Unlike systems with only one component [14j, 
the ratio between the hard and soft cores is critical even if they are comparable to each other since 
different patterns may appear by just changing the value of A. In Fig. 2b we show the structures 
obtained for the same set of parameters as in Fig. 2a except that the value of A has been increased to 
A = 2. The particles now form aggregates in which the particles prefer to align forming intercalated 
stripes of alternating species. Each stripe consists of only one type of particles and try to avoid 
overlap between the coronas of neighboring stripes of the same species, thus, the shoulder width 
act as spacer between neighboring stripes (see inset of Fig. 2b). Increasing even more the range of 
the soft corona, the system develops a completely different strategy to form stable configurations 
giving rise to the appearance of new patterns. This is shown in Fig. 2c where a pattern formed by 
a square array of tetramers develops if A is increased to A = 3. Each square tetramer is formed by 
four close packed particles of the same species and the separation between the tetramers is dictated 



5 



by the condition that the coronas of each particle overlaps at most four coronas of particles of the 
same species, as shown in the inset of Fig. 2c where we have highlighted the coronas that settle 
the spacing between tetramers. Two lines of defects, one near the upper side and the second near 
the right side of the simulation cell are in evidence. They are originated in the periodic boundary 
conditions; effectively, the lattice of tetramers extend over the entire length of the simulation cell 
and join up with their periodic images. Since this length is not commensurable with the periodicity 
of the lattice the defect lines appear. A few vacancies scattered along the lattice that could not 
be annealed out are also present, but this is to be expected for kinetic reasons. For very large 
values of A (A = 10 in Fig. 2e and f), larger hexagonal close-packed aggregates of particles of the 
same species are formed and this aggregates in turn self assemble in square arrays [panel (e)] or 
thick stripes [panel (f)] formed similarly of hexagonal close-packed particles (inset Figs. 2e and 
f). We can think of the aggregates as larger interacting particles with the internal structure only 
changing the effective inter-aggregate interactions [16j, thus the system try to arrange particles so 
that the shape of the cluster becomes as circular as possible. This, in turn, guarantees that the 
underlying structure is close to the energetically most favorable hexagonal lattice [17j. Within the 
lane structure (Fig. 2f) the system follows a similar strategy consisting in first optimize the packing 
inside a stripe, leading to hexagonal particle arrangements inside each lane. Nevertheless, due to 
the increasing complexity of the inner structure of the patterns, a simple energetic explanation in 
terms of overlapping coronas is very difficult. 

In our model, configurational energy is essentially determined by geometry; thus, simple geo- 
metrical considerations can be used to predict the parameters at which a given pattern is expected 
to be adopted at zero temperature. For values of the range of the soft coronas such that A > \pi 
the lattice constant a at T* = is determined by the condition that the soft coronas of the nearest 
neighbors of the same species just touch as shown in the inset of Fig. 2a. Therefore, the lattice 
parameter takes the value a — Actq and the density of the aggregate is Psq^ares — 2/A^. On the 
other hand, if A < \pl then the particle hard cores just touch forming a close packed square array 
whose lattice constant a — v^ctq is determined by the hard-core diameter ctq, and the density of 
the aggregate is Pgquares — The energy per particle of this array is — —46 which is the number 
of overlaps between the corona of a given particle and its first neighbors times the energy of each 
overlap. The square array has the characteristic that no overlap between the coronas of particles 
of the same component takes place. However, if the density is so high that the coronas of particles 
of the same species have not enough space to accommodate without overlapping, then the system 
adopts a different strategy in order to minimize its energy; the formation of lanes provides the 
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energetically best solution. This is achieved (if A > 2 ) when the soft corona of a given particle 
just touch the ones of the second nearest neighbors of the same stripe determining the separation 
between the particles in the same stripe. On the other hand, if A < 2, neighboring particles of 
the same stripe are in direct contact, forming a one-dimensional close-packed arrangement As 
a result, in both cases, the particles form a centered rectangular arrangement as shown in the 
inset of Fig. 2b, with the lattice parameter in the direction of the stripe given by a = Acro/2 if 
A > 2 or a = (Jo if A < 2. The lattice parameter perpendicular to the stripes is given by the 
distance between adjacent stripes of the same species which are separated by a distance b = Actq 
if A > or b = v^ctq if A < \/3. The densities are, correspondingly, p*tripes = 4/A^ if A > 2, 
Pstripes = 2/ A if a/3 < A < 2, and p*tripes = 2/\/3 ~ 1.1547 if A < a/3; the last situation corresponds 
to an hexagonal close packing. Using similar geometrical considerations we obtain the energy per 
particle in each case, u = — 6e if A > ^/3 or u = — 2e if A < \/3. Finally, for the square array of 



tetramers the lattice parameter is a = cro/y/2 (1 + \/2A^ — 1 ] if A > ^/5 and the corresponding 



particle of this array is given hy u = —8e if A > a/5, u = — 4e if 2 < A < a/5, or = if A < 2. 

By comparing to each other the energies of the lattices found previously, it is possible to con- 
struct a phase diagram at zero-temperature which allows to predict the structures adopted for a 
given set of parameters. Regions of the phase diagram where minimum energy configurations at 
r* = appear are shown in Fig. 3 for the three simplest structures obtained in the simulations. 
Here we see, for instance, that stripes are preferred over squares at larger densities and for a larger 
range of the coronas, and that tetramers will be preferred at even larger values of A. The previous 
results, obtained using only geometrical arguments, agree with the structures obtained in Fig.2a, 
b, and c; nonetheless, we do not expect an exact correspondence between the phase diagram at 
T* = with finite temperature simulations since in that case entropy plays an important role. 
This is put in evidence in Fig. 2d which shows the square array adopted in a system with the 
same parameters as in Fig. 2b except that the density is reduced to p* = 0.3. According to the 
phase diagram, for this set of parameters one expects the formation of stripes at T* = 0; however, 
since the square array is less dense, then the particles have more room to vibrate around their 
equilibrium positions which in turn means that the contribution to the free energy F = U — TS 
due to the entropy associated with the thermal motion is more important for the square array than 
for the stripes and therefore the system may adopt the square array. Thus, polymorphic transitions 
are expected under cooling, specially for parameters near the border of the different regions, i.e.. 
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A ^ a/3 and a/5. 

Until now we have focused only in the minimum energy configurations, however, energy consid- 
erations gives no information about the mechanical stability of the structure under the action of 
external agents [2J. The stability of an hexagonal phase (or an fee phase in 3D) is often attributed 
to the optimality of the close-packed lattice [16j. On the other hand, for the square array with 
A > a/2 the touching coronas of the particles belonging to the same species provide mechanical 
stability to the structure; however, this is not the case for shorter ranges of the coronas and the 
structure is unstable with respect to shear. 

Here we are primarily interested in the crystalline configurations that the system may adopt 
at low temperatures. Nonetheless, the behavior of the system at higher temperatures is also very 
interesting since under cooling other structures may preempt the crystalline lattice, as previously 
discussed. In fact, we have found that the original isotropic fluid phase may first suffer a transfor- 
mation to a disordered fluid-like strip phase. This is exhibited in Fig. 4 where snapshots taken at 
two temperatures along the annealing isochore are shown. Starting with an isotropic fluid phase, 
as the temperature is decreased the system first shows the formation of a low-temperature fluid 
phase [panel (a)] consisting of a disordered stripe structure which under further cooling undergoes 
an abrupt transition at a given temperature form the disordered stripe pattern to a crystalline 
square array of tetramers [panel (b)]. For the chosen parameters, the array of tetramers appears 
to coexist with a stripped structure in which part of the configurational disorder typical of high 
temperatures remains. 

The complexity of the patterns grows with increasing asymmetry in the stoichiometry of the 
sample. Two representative configurations for non-equimolar mixtures are shown in Fig. 5. In 
panel (a) we observe a phase separation with a region consisting of a high-density blend of the two 
species and a region consisting of a low-density phase of the majority component. Lanes similar to 
those obtained in the equimolar mixture are observed in the high-density component. In general 
we find that patterns similar to those found for the equimolar mixture appear in the high-density 
phase. The pure component may also form patterns in the surrounding regions. As expected, the 
structures found in this pure phase are also similar to those found in a single-component square- 
shoulder system [l3j,[T3J. A different behavior is observed in Fig. 5b. Here, there is no macroscopic 
phase separation; instead, the system exhibits an ordered hexagonal structure consisting of clusters 
of the minority species corralled by particles of the majority component that form an open network 
structure. It has been proposed that these morphologies can be employed for the fabrication of 
nanostructures [Ij. For example, after selectively removing one of the blocks, the remaining pattern 
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can be transferred into a functional material [9], [20] -[25]. Similar nanoparticle corrals have recently 
been obtained in network-self-assembled monolayer hybrid systems [Ij and are of particular interest 
as surface templates because they contain cavities that can be filled by a variety of guest molecules. 
This kind of open network is particularly flexible since the pore size can be controlled by changing 
the range of the corona and the asymmetry of the mixture. 

It is known that computer simulations risk to be trapped in local energetic minima due to 
the rough and complex energy landscape [17j; thus, in order to circumvent this problem we have 
followed a very accurate simulated annealing procedure since for sufficiently slow cooling down 
of the system one expect that the particles accommodate to its minimum energy configurations 
when temperature vanishes [26j. Nonetheless, alternative approaches could also be of great help to 
predict the complete set of ordered equilibrium structures of the system; one of these approaches 
consists in applying genetic algorithm search strategies to find in a systematic way minimum 
energy configurations [19j. Such algorithms have proven to be a reliable way to predict the ordered 
equilibrium structures for different systems, including monolayers of binary dipolar mixtures [27j. 

In addition to illustrate the complexity of structures that can arise from the simple model 
presented in this paper, the system is possible to be carried out experimentally thanks to the 
recent advances to control colloidal and nanoparticle interactions [8J. From the huge variety of 
different ordered equilibrium structures developed in the system, here we have shown just a few 
of them that were chosen due to their possible practical relevance and because some of them have 
already been observed experimentally [8J. The large versatility of the model indicates that such 
systems can be adapted to a wide variety of technological applications ranging from the directed 
growth of molecular wires or carbon nanotubes [15j to nanolithography and nanoelectricity [13j or 
even as templates for metal deposition. 
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FIG. 1: Description of the model, (a) Two kinds of diblock copolymer, A-B and B'-C. Block B contains 
small numbers of groups (blue) that form hydrogen bonds to complementary groups (red) in block B'. 
(b) Schematic representation of the blend of block copolymer with hydrogen bond interactions forming 
cylindrical micelles (top view), (c) The interaction between a pair of cylindrical micelles is modeled by 
a hard-core soft-corona potential. If the micelles are of the same component then the corona is a square 
shoulder otherwise the corona is a square well. 
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FIG. 2: Spatial arrangements of the system at temperature T* = 0.1 for an equimolar mixture x = 1. (a) 
Square array formed in a system with A = 1.5 and p* = 0.5, (b) stripe phase obtained for A = 2 and p* = 0.5, 
(c) a square array of tetramers appears for A = 2.5 and p* =0.5, (d) square array obtained for A = 2 and 
p* = 0.3. Comparing panels (b) and (d) we see that at the lower density shown the system aggregates in 
square arrays whilst at the larger density the system prefers to form stripes of alternating species. For very 
large values of A, larger aggregates of hexagonal close-packed particles of the same species are formed. These 
aggregates in turn self-assemble in (e) square arrays for A = 10 and p* = 0.3, and (f) thick stripes for A = 10 
and p* = 0.65. 



12 



1.15 
1.0- 

0.8- 

0.6- 



1% 



Squares 

Stripes 

Tetramers 

Others 



0.4- 



0.2- 



Squaj'es StrijDes Tetramers 



0.0- 



2i/2 3I/2 2 5. 



FIG. 3: Zero-temperature phase diagram. This figure shows the different regions of the parameter space 
where different patterns appear. The blue zone represents the region where squares are present; the red 
zone corresponds to the appearance of stripes; finahy, tetrameres appear in the green zone. The yellow zone 
at the upper-right corner of the figure corresponds to other structures not studied in this work. The blue 
(dotted), red (dashed) and green (dash-dot-dot) define the maximum density at which squares, stripes, and 
tetrameres could also appear, respectively, at finite temperatures. 




FIG. 4: Spatial configurations at temperatures T* = 0.5 and T* = 0.1. Snapshots taken at two temperatures 
along the annealing isochore. At high temperatures the system self-assemblies forming fluid-like disordered 
stripes (a), while at low temperatures the system crystallizes in an ordered square array of tetramers (b). 
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FIG. 5: Two representative configurations for asymmetric mixtures {x = 2 and p* =0.5). In (a) we observe 
a phase separation with a region consisting of a high-density blend of the two species and a region consisting 
of a low-density pure phase of only the majority component, when A = 1.5. In (b) we observe the formation 
of colloidal corrals, where the majority component frames clusters of the minority component when A = 2. 



